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ABSTRACT 

We present new radio and optical observations of the colliding- wind system WR 146 
aimed at understanding the nature of the companion to the Wolf-Rayet star and the 
collision of their winds. The radio observations reveal emission from three components: 
the WR stellar wind, the non-thermal wind-wind interaction region and, for the first 
time, the stellar wind of the OB companion. This provides the unique possibility 
of determining the mass-loss rate and terminal wind velocity ratios of the two winds, 
independent of distance. Respectively, these ratios are determined to be 0.20±0.06 and 
0.56±0.17 for the OB-companion star relative to the WR star. A new optical spectrum 
indicates that the system is more luminous than had been believed previously. We 
deduce that the "companion" cannot be a single, low luminosity 08 star as previously 
suggested, but is either a high luminosity 08 star, or possibly an 08+WC binary 
system. 
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1 INTRODUCTION 

Wolf-Rayet stars are surrounded by dense stellar winds giv- 
ing rise to free-free emission extending from IR to radio 
wavelengths. Typically, this emission is characterised by a 
power- law spectrum of the form S v oc v a , with values of 

the spectral index a ~ +0.7 1-0.8, and radio brightness 

temperatures ~ 10 4 K. A small number of WR stars have 
radio emission that exhibits quite different properties: neg- 
ative spectral indices and brightness temperatures ~ 10 6 K 
or higher, properties that are characteristic of non-thermal 
emission. WR 146 is a member of this group, which includes 
WR 125, WR 140 and WR 147. 

The radio emission from WR 146 was first resolved in 
high resolution observations with MERLIN (Dougherty et 
al. 1996, hereafter Paper I). These 5-GHz observations re- 
vealed two components, N5 and S5, separated by ~ 120 milli- 
arcseconds (mas). The flux of S5 was consistent with that 
estimated from extrapolation of the IR-millimetre spectrum 
arising from the free-free emitting envelope around the WR 
star. The brightness temperature of N5 (~ 10 6 K) identi- 
fied the nature of the emission from this component as non- 
thermal. An optical spectrum showed evidence for absorp- 
tion lines at US and H7, which we attributed to an early-type 
companion to the WR star. This led us to hypothesize in Pa- 
per I that the non-thermal emission arose from a population 
of relativistic electrons, accelerated in a wind-wind collision 



region where the wind of the WR star and the companion 
interacted (e.g. Eichler & Usov 1993). To be consistent with 
such a model, we suggested the companion lay at the same 
position angle as N5 from S5, but slightly further away from 
the WR star. 

The presence of a companion was confirmed in opti- 
cal imaging with the Hubble Space Telescope (HST) by 
Niemela et al. (1998). They observed two stars, WR146A 
and B (hereafter So and No respectively), at the same posi- 
tion angle as the radio sources but separated by ~ 168 mas. 
Under the assumption that the southern sources in both the 
HST and MERLIN images are coincident, these observa- 
tions place the non-thermal source between the two stellar 
images, strongly supporting wind-wind collision as the ori- 
gin of the non-thermal emission. Its position relative to the 
two stellar components (~ 120 mas from So and ~ 48 mas 
from No) is where the dynamical pressure of the two stel- 
lar winds is balanced. This indicates that the momentum of 
So's stellar wind is ~ 0.1 times that of No- With the wind 
velocity of the WC star in WR 146 (~ 2900 km s~\ Eenens 
& Williams 1994) being greater than that of a typical OB 
star, and the expectation that the mass-loss rate of a WR 
star would be greater than that of an OB star, this strongly 
supports the identification of So with the WC6 star in Pa- 
per I and identifies No with an OB companion having a 
lesser wind momentum. The momentum ratio lead Niemela 
et al. to infer that the wind momentum of the companion 
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was more appropriate to a star of early O or Of type than a 
late O-type main-sequence star. Taking account of the pho- 
tometry, they suggested an 06-05 V-III spectral type for 
the companion. 

The photometry by Niemela et al. showed the two stars 
to be equally bright in B but that So was redder than No 
in (B-V) and (U-B). The (near-zero) magnitude difference 
between the two blue images falls between two very different 
estimates of the WR:0 light ratio deduced from spectra of 
WR 146 which included both stars. 

On the one hand, in Paper I, we measured equivalent 
widths of H<5 and H7 in our blue spectrum of WR 146 to 
be W\ ~0.9A. Comparison of these with those (~ 2.5A) 
typical of mid-to-late O stars in the Walborn & Fitzpatrick 
(1990) atlas indicated that the O-star spectrum was diluted 
and suggested that the continuum of the WC star was twice 
as bright as that of the O companion in the blue. On the 
other hand, Willis et al. (1997) used a spectrum of lower 
resolution but longer wavelength coverage than presented in 
Paper I to determine the WR:0 light ratio from comparison 
of the equivalent widths of the emission lines with those 
of two single WC6 stars. They found that the O star was 
brighter than the WC star, reporting a continuum light ratio 
(WR:0) of 1:(2±1). 

The discrepancy in the light ratios of the WC6 and O 
components derived from the dilution of O star absorption 
lines and WC6 star emission lines is significant and probably 
too great to arise purely from the observational uncertain- 
ties. Leaving aside the question of whether one of the stars 
varied between the two observations, the principal uncer- 
tainties in two spectroscopic light-ratio determinations come 
from whether the intrinsic line strengths of the components 
are indeed equal to those of the comparisons adopted. The 
discrepancy could be reduced if the equivalent widths of H<5 
and H7 in the O star were weaker than the 2.5A adopted or if 
the emission lines in the WC6 component were weaker than 
those of other WC6 stars. Both possibilities are plausible. 

The strengths of the emission lines in WR 146 may well 
be atypical for its WC6 type: both Eenens & Williams (1992) 
and Willis et al. found anomalously low C/He abundances 
for this star. On the other hand, the strengths of H7 and H<5 
in O-type stars do depend on spectral type and luminosity 
class, which are not known directly for the companion in 
WR 146. We note that the "08.5V" adopted by Willis et al. 
comes not from the star's spectral lines but from its luminos- 
ity inferred from that of the WC6 star and the continuum 
light ratio. We therefore re-observed the blue spectrum of 
WR 146 with the William Herschel Telescope in an attempt 
to determine the spectral type and, hopefully, luminosity of 
the early-type companion to the WR star. These data, to- 
gether with the photometry by Niemela et al., will give us a 
better idea of the intrinsic properties of the companion for 
comparison with the radio observations. 

We also re-observed WR 146 at radio wavelengths, ex- 
tending the frequency coverage of Paper I, to search for the 
wind flux from the companion, and study the characteristics 
of the two previously observed radio components with the 
aim of furthering our understanding of the colliding wind 
phenomenon. 
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Figure 1. a) 5-GHz MERLIN observations obtained on 1992 De- 
cember 26 The la rms uncertainty is 250/iJy, with contour levels 
at -3, 3, 6, 9, 12, 24, 48a. The synthesized beam is 53 X 53 mas 2 . 



2 RADIO WAVELENGTH STUDY 
2.1 The observations 

Multi-frequency observations of WR146 at 1.4, 5, 8.4, and 
22 GHz were taken on 1996 October 26 using the A- 
configuration of the NRAO Very Large Array (VLA). Ob- 
servations of the nearby radio-bright quasar 2005+403 were 
interleaved with the observations of the target source for 
phase-referenced calibration of the antennae gains. The ab- 
solute flux scale was determined by observation of 3C48, 
assumed to have fluxes of 15.970, 5.516, 3.226 and 1.174 Jy 
at 1.4, 5, 8.4, and 22 GHz respectively. 

In addition, archival data from MERLIN at 5 GHz have 
been analysed. A 4-hour on-source observation at 5 GHz 
was obtained with MERLIN on 1992 December 26, along 
with frequent observations of the phase calibration source 
2005+403. Following the standard procedure for calibrating 
MERLIN 5-GHz data, the flux scale was established by ob- 
servation of 3C286 with an assumed total flux of 7.309Jy 
(7.020 Jy on the MK2-Tabley baseline) and the unresolved 
source 0552+398, which had a derived "bootstrap" flux of 
5.786 Jy. The bootstrap fluxes for 2005+403 are given in 
Table [if 

It is possible to check the accuracy of the derived flux 
scales in the 5- and 8.4-GHz observations. Comparison of the 
"bootstrap" fluxes of 2005+403 at these frequencies to those 
given in the University of Michigan Radio Astronomy Ob- 
servatory (UMRAO) flux database shows that, where closely 
contemporaneous observations are available, the agreement 
with UMRAO fluxes at these frequencies is within 1%. Given 
the excellent agreement at these two frequencies, it is rea- 
sonable to assume that the flux scales at the other observed 
frequencies are good. 

Aside from the initial amplitude calibration of the 
MERLIN observations, the calibration and subsequent 
imaging of the data were carried out using the NRAO AIPS 
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Figure 2. VLA A-configuration observation at 22 GHz obtained 
on 1996 Oct 26. Two components are clearly resolved. The la 
rms uncertainty is 0.5mjy, with contour levels at -3,3,6,9,12<r. 
The synthesized beam is 100 X 75 mas 2 at a position angle of 
-2°. 



Table 1. Radiometry of WR146 



VLA observations 


Epoch 


V 


N 22 


S22 


2005+403 




(MHz) 


(mjy) 


(mjy) 


(Jy) 


96/10/26 


22460 


10.4 ± 1.0 


7.0 ± 1.3 


2.136 




8435 


29.8 ±0.8 


2.956 




4885 


37.6 ±1.0 


3.227 




1465 


78.4 ±0.2 


2.596 


MERLIN observations 


Epoch 


V 


N 5 


s 5 


2005+403 




(MHz) 


(mjy) 


(mjy) 


(Jy) 


92/12/26 


4885 


31.4 ±0.4 


2.5 ±0.3 


2.923 


95/04/29 


4885 


28.5 ±0.3 


1.6 ±0.3 


2.865 



Notes: The parameters for WR146 at epoch 95/04/29 have been 
re-determined from the original data. The flux values are slightly 
higher, although consistent within the uncertainties, than the 
values originally quoted in Paper I. The quoted uncertainties in 
the fluxes are the formal-fit errors. The uncertainty in the 
absolute flux scale (~ a few %) is not included and, for at least 
the higher fluxes, will be an important source of uncertainty. 



software package. The complex antennae gains were ini- 
tially derived for 2005+403 and then interpolated (phase- 
referenced) to the observations of WR 146. In addition, sev- 
eral iterations of phase-only self-calibration were used to re- 
fine the antennae gains during the MERLIN observations of 
WR 146 to improve the dynamic range of the final synthe- 
sized images. 

At 1.5 and 4.9 GHz the VLA data reveal an unresolved 
source, whereas at 8.4 GHz the radio emission is marginally 




wrl46s (23GHz) 




Baseline distance (MA) 

Figure 3. The 22-GHz real visibility distributions for N22 (top) 
and S22 (bottom). The data were binned into nine annuli in the 
uv-plane, and then averaged. The averages and their uncertain- 
ties are shown versus baseline distance. The imaginary visibility 
distributions are consistent with zero. The solid lines represent 
weighted circular Gaussian model fits to these data. The FWHM 
of the fits are 48 ± 11 and 78 ± 19 mas respectively. It is clear 
from these figures that these data do not merit more sophisticated 
modelling. 



resolved. However, WR 146 is well resolved into a double ra- 
dio source at 22 GHz. The final synthesized 22-GHz image 
is shown in Fig. |^. We identify the northern and southern 
components as N 22 and S 22 . The visibility data suggest both 
these components are resolved (Figure pf) . Gaussian model 
fits to the visibilities give the diameter of the radio emitting 
regions as 48 ± 11 and 78 ± 19 mas respectively (see Paper 
I for a description of the method). The large relative un- 
certainty in these sizes is a direct consequence of the high 
rms uncertainty per visibility at 22 GHz. The fluxes from 
the VLA observations were determined by Gaussian source 
fitting using the AIPS routine JMFIT. The measured fluxes 
are given in Table [[]. 

The final synthesized MERLIN image at 5 GHz is shown 
in Fig [l|. The image reveals two components of emission, 
identified as N5 and S5. In the MERLIN observation, the 
flux and size of component N5 were determined by Gaussian 
model fitting of the visibility data. The derived diameter for 
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Table 2. J2000 Positions of emission components 



20 33 59.10 59.09 59.08 

RIGHT ASCENSION (B1950) 

Figure 4. Overlay of the VLA 22-GHz image (thin contours) 
with the MERLIN 5-GHz image from 1992 (thick contours), as- 
suming S22 and S5 are coincident. Alignment was attained using 
the brightest point in the two southern components. Clearly, the 
northern component at 22 GHz is displaced to the north of the 
non-thermal region observed at 5 GHz. We identify N22 as the 
stellar wind emission originating from the OB companion to the 
WC6 star. 



N5 is 42 ± 1 mas. Gaussian fitting to the image data were 
used for the source parameters of S given in Table 1. 

Qualitatively, the 1992 MERLIN image is consistent 
with the 5-GHz observation from 1995 April 29 (Paper I) 
though there is some evidence that N5 may have decreased 
in flux by ~ 3 mjy between 1992 and 1995. The flux of S 5 
appears to have decreased between the two epochs. However, 
the difference is only 2a, much less than the 5a threshold 
typically adopted as evidence for variation. Therefore, we 
adopt the mean value of 2.0 ± 0.2 mjy for the flux of S5. 
Variations in the total 5-GHz emission from WR 146 have 
been observed using the WSRT over the last decade that 
are attributed to variations in the non-thermal component 
(Setia-Gunawan et al. 2000). Both our MERLIN 5-GHz to- 
tal flux values are consistent with the WSRT observations. 
Variability could account for the slightly higher 5-GHz flux 
at the epoch of the VLA observations. 

The positions of the components observed at 5 GHz 
and 22 GHz are quoted in Table ^[ The absolute position of 
the northern component in the various observations was de- 
duced from the phase-reference only calibrated images. Ab- 
solute position information is lost during the self-calibration 
process, though the relative position is preserved. These 
were deduced from the images shown in Fig. ^ and Fig. [| 



Epoch 


Source 


RA 


Dec 


A 


PA 






2Qh 35m 


41° 22' 


(mas) 


(°) 


fir* / -1 / 

92/12/26 


N(5) 


47 .0905 


44 .832 








S(5) 


47.0885 


44.713 


122 ±4 


14 ±3 


95/04/29 


N(5) 


47.0896 


44.811 








S(5) 


47.0854 


44.705 


123 ±4 


31 ±3 


96/10/26 


N(22) 


47.0997 


44.754 








S(22) 


47.0959 


44.603 


162 ± 8 


22 ±4 



Notes: the data from Paper I have been re-analysed giving the 
new absolute positions quoted here. The accuracy of the 
absolute positions is ~ ±10 mas, due to the uncertainty in the 
optical position of the phase-reference calibrator 2005 + 403. 
However, the relative position error is largely determined by our 
ability to fit a Gaussian function to the emission from 
WR146(S). 



2.2 Identifying the radio components 

The nominal resolution of A-configuration at 22-GHz obser- 
vations is very similar to that of MERLIN at 5 GHz, allowing 
direct comparison of these data. At 22 GHz, the VLA obser- 
vations reveal two components N22 and S22, very similar in 
appearance to the MERLIN 5-GHz data. However, the rela- 
tive positions of the two 22-GHz components show them to 
be significantly further apart than those from the MERLIN 
5-GHz data, leading to the question of whether the 22-GHz 
sources are the same as those observed at 5 GHz. 

The brightness temperatures of N22, S22 and S5 are 
1.6 x 10 4 K, ~ 5000 K and ~ 4000 K respectively, all consis- 
tent with an origin in a photo-ionized circumstellar envelope 
where the photo-ionization equilibrium temperature is typi- 
cally ~ 10 4 K. In contrast, this is about two orders of mag- 
nitude lower than the brightness temperature of 1.3 x 10 6 K 
for N5, which clearly indicates a non-thermal origin. 

The separation of S22 and N22 (162 ± 8 at a position 
angle of 22 ± 4°) is very close to that (168 ± 31 mas at posi- 
tion angle 21 ± 4°) of the optical components (No and So) 
observed with the HST. There is no evidence of intrinsic 
proper motion in the source at 22 GHz in recent (1999) ob- 
servations of the source (A. Fink, private communication). 
The coincidence of the optical and 22-GHz positions, the 
lack of evidence of intrinsic proper motion, and the thermal 
nature of the emission from the two 22-GHz components 
leads us to conclude that the 22-GHz components are the 
stellar winds associated with the two stellar components im- 
aged with the HST by Niemela et al. (1988). This is the first 
time that the stellar wind of the companion in a WR binary 
has been spatially resolved. 

In Paper I we concluded that S5 was the thermal emis- 
sion from the WR-star wind. Indeed, if we take the mean 
flux value for S5 from our two epochs of MERLIN obser- 
vations, the spectral index between 22 and 5 GHz for the 
southern component is +0.82 ± 0.14, consistent with those 
observed in the stellar winds of other WC-type WR stars 
(e.g. Williams 1996). We estimate the diameter of S5 to be 
~ 190 mas from the diameter of the S22, since angular size 
0.8. This is approximately four times larger 



© 0000 RAS, MNRAS 000, 000-000 



WR 146 - Observing the OB-type companion 5 



than that given in Paper I. However, we feel the diameter 
presented here is a more reliable estimate. 

Additionally, in Paper 1 we identified N5 as non-thermal 
emission from a wind- wind collision region. If we overlay the 
22-GHz and 5-GHz data assuming two southern components 
originate in the WR star wind (see Figure 4) we can see that 
the relative separation of the components supports such a 
model, where the wind-collision region must fall between the 
two stars. 



2.3 Wind parameters for both stellar components 

For a steady state, smooth, fully ionized stellar wind having 
a radial r~ 2 ion density distribution, the free-free radio flux 
S u is related to the stellar-wind density i.e. M/v by 



S v = 2.32 x 10 4 (7g„,T e ^) 1 



( 



d~ 



mjy, 



(1) 



where v is the frequency in Hz, d is the distance in kpc, 
M is the mass-loss rate in M© y , v the terminal velocity 
of the wind in km s _1 ; 7, Z, fi and g v ,r e are respectively, 
the number of electrons per ion, the mean charge per ion, 
the mean atomic weight and the Gaunt free-free factor, a 
function of electron temperature and frequency. It follows 
from that 



M 



53. 



(2) 



Thus, the measurement of radio emission from the two stel- 
lar wind components in the WR 146 system leads to the den- 
sity ratio of the two winds. If the stellar winds are clumped 
rather than smooth, we can allow for this by replacing M in 
the above with (M/^/f), where / is the wind filling factor, 
assumed to be constant over the radio-emitting region. 

In the case of the collision of two spherical winds at 
terminal velocity, the contact discontinuity between the two 
winds intersects the line of centres between the stars at the 
point of momentum balance. If the projected separation of 
the WR and OB star is Dcosi, and that of the OB com- 
panion from the non-thermal region is tob cos i, we can then 
write 



roti/D = 



1 + 7?" 



(3) 



where 77 is the ratio of WR and OB-companion wind mo- 
menta, (Mj))ob/(M«)wr- The wind-collision geometry and 
equation |^ would also be affected by clumping, but this is 
beyond the scope of the present study. 

By defining the ratio of the OB star and WR-star wind 
densities, (M/u)ob/(M/u)wr, as £, then 



Mqb 
M WR 



& 



VQB _ 



(4) 



Knowing Mwr and uwr we can now uniquely determine 
Mob and i>ob that satisfy simultaneously both £ and 77. An 
attractive property of these ratios is that they are indepen- 
dent of distance, typically a very uncertain parameter for 
WR stars and particularly for WR 146. 

For the OB companion we will assume that hydro- 
gen and helium are, respectively, singly and doubly ionized, 



leading to values 7 = 1.1, Z = 1.15, /i = 1.34 (Lamers 

6 Leitherer 1993) and T e ~ 10 4 K. Taking the values of 

7 = 1.15, Z = 1.2, fj, = 5.29, T e = 8 000 K from Willis 
et al. (1997) for the WR star, and that the fluxes of N 22 
and S22 are those of the OB star and the WC star respec- 
tively, implies that £ = 0.36 ± 0.19. From Table 2, D is 
162 ± 8 mas and tob is 40 ± 9 mas, so equation ^ gives 
77 = 0.11±0.03. Thus, the ratios of mass-loss rates and wind 
velocities are 0.20 ± 0.06 and 0.56 ± 0.17 respectively. The 
largest uncertainty in these ratios arises from £ . Though this 
requires knowledge of the relative metallicity and ionization 
structure within the two winds, the values of Z and 7 are 
closely the same and effectively cancel. On the other hand, 
the mean molecular weights of the two winds are quite dif- 
ferent and contributes a factor ~ 2 to the ratio. The relative 
uncertainty in the observed fluxes provides the bulk of the 
uncertainty in £. 

From the deduced velocity ratio and the terminal wind 
velocity measured for the WC6 star (2900 km s _1 , Ee- 
nens & Williams 1994), we derive a terminal velocity of 
the OB wind of 1600 ± 480 km s" 1 . This is consistent with 
those determined for late O-type stars of any luminosity 
class (Prinja, Barlow & Howarth 1990). However, the de- 
duced ratio (0.19) of mass-loss rates implies a very high 
mass-loss rate for the OB star, irrespective of the distance 
to WR 146. For example, if the WC6 star has a mass-loss 
rate comparable to the average rate derived from radio or 
from "Standard Model" analyses of optical emission lines 
(4 x 10~ 5 M y" 1 , Willis 1999), that of the OB star would 
be ~ 8 x 10 -6 M0y _1 . This is comparable to that deter- 
mined for the 08f supergiant HD 151804 (Bieging, Abbott 
& Churchwell 1989, Crowther & Bohannan 1997) and al- 
most 50 x greater than the mass- loss rates typical of main- 
sequence 08 stars (~ 1.7 x 10~ 7 M y" 1 , Howarth & Prinja 
1989). This result still holds if the true mass- loss rates of 
WR stars are lower than the value cited owing to clumping 
since it is really the ratio of (M/y/f)s that has been de- 
termined. Only if the OB-star wind was clumped would its 
mass-loss rate be reduced by a factor of y/f . 

Throughout the paper we have identified the southern 
and northern sources as the WR and OB star respectively. 
If these identifications were switched i.e. WR star is to the 
north, and the OB star to the south, then the deduced wind 
velocity of the OB star is greater than 18 000 km s _1 and 
the mass-loss rate is 1.4 times that of the WR star. Clearly 
these parameters are unacceptable given our knowledge of 
OB-star winds and we dismiss such a source identification 
as untenable. 

The mass-loss rates determined from the measured free- 
free wind fluxes depend heavily on the adopted distance to 
WR 146. Therefore, we first re-examine the luminosity of the 
system using new optical observations. 



3 OPTICAL SPECTROSCOPY 
3.1 The observations 

The spectrum of WR 146 was taken with the blue arm of 
the ISIS spectrograph on the William Herschel Telescope 
(WHT) on 1999 June 19. As in Paper I, we observed in the 
blue, where the relative contribution of the OB companion 
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Figure 5. Rectified blue spectrum of WR146 with absorption features marked. Above it is the spectrum of another WC6 star, WR15. 
This has been artificially diluted so as to match the emission-line spectrum of WR 146 as an aid to determination of the continuum 
in WR 146 for measurement of the equivalent widths of the absorption features. Some of the diluted WR emission features are marked 
above the WR 15 spectrum. 



would be greatest, and avoided the strong A4650A emission 
feature from the WR star. The R1200B grating and 1.1 arc- 
sec slit gave a resolution of 0.9 A. Four integrations, each of 
750s, were taken of the star, interleaved with observations 
of an argon discharge tube for wavelength calibration. The 
detector used was a 4200 x 2048 EEV42 device. 

We present the new spectrum of WR 146 in Fig. ^ and 
list absorption lines in Table |^. The spectrum has been 
smoothed giving a resolution ~ 1.0 A (FWHM), as mea- 
sured from the widths of the narrow interstellar CH, CH + 
and Cal lines. The FWHM of the A3933A Call line is 1.6 
A, the greater width coming from optical depth effects (cf. 
Jenniskens & Desert 1994). The well developed interstellar 
spectrum is consistent with the heavy reddening suffered by 
WR 146. 

The broad undulations in the continuum are mostly 
emission lines from the WC6 star, diluted by the light from 
the companion. The emission lines in WR 146 are unusually 
broad for its WC6 subtype (cf. Eenens & Williams 1994, 
Willis et al. 1997), making it easier to measure the absorp- 
tion spectrum of the companion. The absorption line widths 
are typically ~ 5 A FWHM, comparable for example to those 
of the 09.5II star 8 Ori (cf. Voels et al. 1989), which has a 
relatively high v sin i (144 kms" 1 , Howarth et al. 1997). 

There is one absorption feature from the WC6 star; 
the P-Cygni absorption component of the A3889A He I line, 
blue-shifted by 2630 kms" 1 to 3854.5A. This line arises 
from the same metastable level as the A10830A line from 



which Eenens & Williams (1994) derived a terminal wind 
velocity (v<x>) of 2900 kms -1 by profile fitting and, like 
the A10830A line, is expected to form in the outer regions 
of the WC6 wind. These results are consistent with the 
Voo = 2700 kms -1 measured by Willis et al. from the [Nelll] 
A15.5-/xm fine-structure line, also formed in the outer wind. 

The other absorption lines (Table |^) are formed in the 
OB-companion star. Each of the hydrogen lines is blended 
with the nearby member of the Pickering series of He II but 
the relative weakness of the unblended odd-numbered Pick- 
ering line at 4200A and the closer proximity of the hydrogen- 
line wavelengths to the observed wavelengths suggests that 
the Hell contribution is relatively small. Helium I is repre- 
sented by the AA 4472 and 4026A lines; the latter is blended 
with He n and is therefore not very useful for diagnostic pur- 
poses. Other He I lines in the Walborn & Fitzpatrick (1990) 
atlas of hot-star spectra are not seen in our spectrum. The 
relative strength of the A 4070A C ill triplet is a little surpris- 
ing but the continuum here is uncertain as it falls between 
two emission features from the WC6 spectrum. 

In order to establish the continuum in the regions of 
each of the absorption lines so that we could measure their 
equivalent widths, we compared our rectified spectrum with 
that of another WC6 star, WR15 (HD 79573), artificially 
diluted so as to match the spectrum of WR 146. The spec- 
trum of WR 15 was observed with the 1.9m telescope at the 
South African Astronomical Observatory and will be dis- 
cussed elsewhere. As a by-product of this matching, we es- 
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Table 3. Measured wavelengths (A) and identifications of absorp- 
tion features in the spectrum of WR 146. The "A lab" values for 
the interstellar features come from Jenniskens & Desert (1994). 



A obs A lab Identification (and equivalent width, Wj) 



3835.4 


3835. 


.4 




3854.5 


3888. 


6 


He I from WC6 star at -2630 km s" 1 


3889.6 


3889. 





HC (?+ 3888.6 He I 2s-3p) 


3933.4 


3933. 


A 


Interstellar Can 


3968.3 


3968. 


5 


Interstellar Call 


3970.4 


3970. 


1 


He 


4025.9 


4026. 


2 


He I 2p-5d + 4025.6 Hell 


4071.1 


4070. 


3 


C in 4f-5g 


4089.6 


4088. 


9 


Si iv W A = 0.06±0.03A 


4102.6 


4101. 


8 


Ho (+4100.0 Hell) 


4200.8 


4199. 


8 


Hen; W A = 0.30±0.03A 


4226.6 


4226. 


.7 


Interstellar Cal 


4232.3 


4232. 


.5 


Interstellar CH+ 


4299.9 


4300. 


3 


Interstellar CH 


4340.9 


4340. 


5 


H 7 (?+4338.7 Hen); W A = 0.70±0.05A 


4430 


4428. 


9 


DIB "A4430" 


4472.1 


4471. 


5 


He I 2p-4d; W A = 0.50±0.04A 


4501.4 


4501. 


.8 


DIB 



timated the dilution of the WC6 spectrum in WR 146, find- 
ing OB:WR « 3, consistent with the ratio, OB:WR =2±1, 
found by Willis et al. (1997) from stronger emission lines be- 
tween 4660 A and 6560 A. We are chary of putting too much 
weight on this result as it depends heavily on the rectifica- 
tion process and the presumption that the WC6 emission 
lines in WR 146 have the same strengths as those in WR 15. 
Instead, we concentrate on the absorption lines, to be dis- 
cussed below. 

3.2 Spectral type and luminosity of the OB 
companion 

To investigate the nature of the OB companion, we begin by 
using the helium-line ratios to estimate its spectral type. Our 
optical spectrum includes contributions from both No and 
So but the ratios of the absorption lines in the companion 
will be equal to those measured from our combined-light 
spectrum, and independent of the dilution by the WR star, 
provided that the dilution does not change significantly over 
the wavelength range of the spectrum. 

Formal O-star spectra classification is based on the ra- 
tio of the A454lA Hen and A4472A Hei lines. Since the 
former was not available, we used the A4200A Hell line 
from the same series. We formed an empirical calibration 
of A4472A/A4200A ratios against spectral type using the 
line strengths measured for a large (~ 100) sample of O 
stars by Conti & Alschuler (1971) and Conti (1973). Using 
this calibration, the ratio of our measured equivalent widths 
(Table ^) indicates a type of 08, with an uncertainty of half 
a subclass. The type is close to that inferred by Willis et al. 
from their OB:WR light ratio and an average luminosity for 
the WC6 star; but this agreement is fortuitous. 

The luminosity class criterion adopted by Conti & 
Alschuler for the middle and late-type O stars was the ra- 
tio of the Si IV A4089A and He I A4143A lines. Unfortu- 
nately, the silicon line is too weak (Wa ~ 0.06A) to mea- 
sure with confidence and the A4143A line is not seen at 



all. Another significant He I line that is apparently miss- 
ing is the A4388A singlet line. The strength of this relative 
to that of the A4472 A He I line ( "singlet to triplet ratio" ) 
has long been known to be sensitive to luminosity owing 
to the relative overpopulation of the He I 2 3 P° state in ex- 
tended atmospheres (e.g. Voels et al. 1989 and references 
therein). Using the equivalent- width measurements by Conti 
& Alschuler (1971) and Conti (1973), we examined the ra- 
tios, W'= W A (A4388)/W A (A4472), formed from their ob- 
servations of the 07.5, 08 and 08.5 stars in their sample. 
This ratio was found to be a strong diagnostic for Of stars. 
All the 07.5-8.5 stars having W' < 0.25 are Of stars, most 
of them supergiants. Measurement of the A4472A line and 
examination of our spectrum in the region of the A4388A 
line indicates that W' <C 0.2, strongly suggesting that the 
08 star in WR146 is an Of star, and possibly a supergiant. 
Unfortunately, the classical Of diagnostic emission features 
Cm AA4630-41A and Hen A4686A fall within the strong 
A 4650A Cm-iv/Heii emission feature from the WC6 star, 
which extends from 4620A to 4690 A in WR 146 (Dessart et 
al. 2000) . This increases the dilution of the OB spectrum to 
~ 1:4 in this critical region and we cannot make a formal Of 
classification. 

Further circumstantial support for a high luminosity for 
the OB star comes from comparison of the A4472A He I to 
H7 ratio with those derived from the Conti & Alschuler and 
Conti (1973) datasets. On the other hand, we do not ob- 
serve emission lines corresponding to the 'unidentified' fea- 
tures observed in 06-09.7 supergiants, including the 08If 
star HD 151804 (Crowther & Bohannan 1997), at 4486Aand 
4504 A. Also, the weakness of the Siiv A4089A line is sur- 
prising given the other evidence for high luminosity, but not 
unique: the Si IV line is barely visible in the spectrum of the 
07Iaf star Sanduleak 80 (Walborn & Fitzpatrick 1990). Al- 
though we cannot assign a luminosity class, the spectrum 
points to an extended atmosphere and high luminosity, but 
not as high as that of HD 151804. 

3.3 Spectral line dilution and the HST 
photometry 

With a better idea of the nature of the 08 star, can we recon- 
cile the measured absorption-line strengths in our composite 
spectrum with the observations by Niemela et al. (1998)? We 
follow Paper I and Niemela et al. by assigning the northern 
optical image (No) to the 08 star on the basis of its closer 
proximity to the non-thermal radio source and the assump- 
tion that its stellar wind has a lower momentum than that 
of the WC6 star. If we assume that the spectrum of So is 
that of a WC6 star, we can estimate the intrinsic equiv- 
alent widths in the 08 spectrum from those observed in 
our combined-light spectrum and Niemela et al.'s photome- 
try of So and No- We find 1.4A for the intrinsic equivalent 
width of H7, which is closer to the average of those measured 
by Conti (1973) for 08f stars (1.3±0.3A) than for 'non-Of 
O8III-V stars (2.5±0.2A), another indication of high lumi- 
nosity properties for the OB star. Our earlier view (Paper 
I) that the WC6 star was brighter than the OB star was 
based on the erroneous assumption of intrinsic absorption- 
line strengths appropriate for a main-sequence OB star. 

The observation of approximately equal B magnitudes 
for No and So is not consistent with the light ratios derived 
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from dilution of the emission-line spectrum of the WC6 star 
by Willis et al. (1997), and supported by our estimate in 
Section 2.1. This could be resolved also if the emission lines 
in the WC6 component were 30-40% weaker than those typ- 
ical of WC6 stars. This is plausible; as already noted, the 
star appears to have an abnormally low C/He abundance ra- 
tio and the dilution analysis by Willis et al. depends heavily 
on C ill and C IV lines. Therefore, the very different light 
ratios originally derived from combined-light spectra in Pa- 
per I and by Willis et al. can be reconciled with a model in 
which WR 146 comprises two stars of approximately equal 
brightness in the blue and relatively weak lines, those of No 
being consistent with other indicators of its high luminosity. 



3.4 The nature of the northern component 

The H/y strength (together with HST photometry), the He I 
singlet-to-triplet line ratio, and the very high mass-loss rate 
all point to an extended atmosphere for No . If it is single, it 
cannot be a main-sequence star but must be luminous. Al- 
ternatively, it may itself be an unresolved binary comprising 
an 08 star with a third component. 

If No was a single, luminous star, we can estimate its 
absolute magnitude from those of 08I/f and 08III stars 
(M v = -6.5 and -5.5, Conti 1988; Vacca, Garmany & 
Shull 1996). For example, combination of an intermediate 
luminosity, Mv = —6.0, which might be appropriate given 
the apparent absence of the 'unidentifed' 4486A and 4504A 
emission lines, with the photometry by Niemela et al. would 
give Mv = —6.2 for So and a combined Mv = —6.9 for the 
WR 146 system. The implied luminosity of the WC6 star, 
Mv = —6.2 or M v ~ —6.0 on the narrow-band ubv system 
used for WR stars, is significantly greater than those of other 
WC6 stars (e.g. mean M v = —4.1 ± 0.4, van der Hucht, in 
preparation), but this may be another manifestation of the 
anomalous nature of the WC6 star in WR 146. 

Adopting a reddening of Eb-v = 2.7, consistent with 
Paper I and Willis et al. (1997), the combined Mv implies a 
distance modulus m — M — 11.2, greater than our distance 
in Paper I (m — M — 10.4) but consistent with member- 
ship of the Cyg OB2 Association (m — M = 11.2, Massey 
& Thompson 1991; Torres-Dodgen, Tapia & Carroll 1991). 
Given the angular distance (~ 0.5°) of WR 146 from the cen- 
tre of the association, stronger evidence is needed to confirm 
membership. 

At this distance (1.7 kpc), the mass-loss rates derived 
for the WC6 and 08 stars are very high: 1.3 x 10~ 4 M y _1 
and 2.5 x 10~ 5 Moy -1 respectively. The latter is three times 
that of HD 151804, whereas the luminosity of No appears 
to be lower. Such a system would have stronger emission 
lines that those in HD 151804, contrary to observations (e.g. 
4486A and 4504A or Ha in Dessart et al. 2000, fig 4 com- 
pared with Crowther & Bohannan 1997), so the single high- 
luminosity model for No is problematic. 

If we consider No to be an unresolved binary in which 
some of the mass-loss is contributed by a third component, 
the latter is unlikely to be another O-type star. Since M oc 
L 1 ' 7 for O stars (Howarth & Prinja 1987), we expect the 
mass-loss rate of an O+O system to be lower than that of 
a single star having the same total luminosity. Also, the 
intrinsic strengths of the H7 lines from two lower luminosity 



stars are likely to be greater, in conflict with our combined- 
light spectrum and the HST photometry. 

A model wherein the third component is another WC 
star suffers neither of these problems. Its spectral subtype 
would have to be similar to that of So, or it would have to 
be fainter, not to have been detected in the combined-light 
spectra. Its presence could account for some of the anoma- 
lous properties (e.g. low C/He) of the WC6 star reported 
in the studies referred to. It could also modulate the wind 
of the 08 star flowing towards the wind interaction region, 
providing a mechanism for the 3.38-y variability in the non- 
thermal emission reported by Setia Gunawan et al. (2000). 
To examine this model, separate spectra of No and So , e.g. 
with an adaptive optics system, are needed to see whether 
No has a WC companion (or Of emission lines) and to de- 
termine its spectrum without dilution. 



4 CONCLUSIONS 

From high-resolution observations with the VLA at 22 GHz 
and with MERLIN at 5 GHz we have observed all three 
components of the WR 146 system: the OB and WC6 stellar 
winds and the non-thermal source where they collide. The 
source geometry and ratio of stellar wind fluxes allow us to 
determine the ratios of mass-loss rates and wind velocities 
independent of distance to the system. From these ratios 
and the observations of the WC6 star, we derive the wind 
velocity of the OB star to be 1600 ± 480 kms -1 and its 
mass-loss rate to be one quarter that of the WC6 star. If 
the WC6 star has an "average" WR-star mass-loss rate of 
~4x 10~ 5 M y -1 , that of the OB star would be ~ 8 x 
10~ 6 M Q y _1 , suggesting a very high luminosity object. 

Support for high luminosity comes from the optical 
spectrum of WR 146, which includes both stars. This shows 
absorption-line ratios formed in the OB companion suggest- 
ing it to be a high-luminosity 08 star. If it is a single star, 
the inferred luminosity places WR 146 at the distance of the 
Cyg OB2 association. This gives an anomalously high lu- 
minosity for the WC6 star but, given its other anomalies, 
does not rule out this distance. However, the mass-loss rates 
determined using this distance (1.3 x 10~ 4 Mqy^ 1 for the 
WC6 star and 2.6 x 10~ 5 M Q y _1 for the 08 star) are awk- 
wardly high and the latter is probably inconsistent with the 
spectroscopy. Many, if not all, of the observations could be 
explained if the companion was itself a binary comprising 
an 08 star and another WC star. This needs to be tested by 
separate spectra of the two visual components of WR 146. 
The presence of an unresolved companion to the 08 star 
could modulate its wind so as to cause the 3.38-y variability 
in the non-thermal emission reported by Setia Gunawan et 
al. (2000). 

We may not have reached a firm conclusion as to the 
nature of the stellar companion(s) to the WC star from the 
new radio and optical observations. However, the ability to 
study all three radio components of the WR 146 system sep- 
arately will ensure that this system becomes an archetype 
for studying the wind-collision phenomenon. 
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